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ABSTRACT 

This  work  is  concerned  with  the  optimization  of  an  ac¬ 
tive  sonar  or  radar  transmit  signal,  to  maximize  the  prob¬ 
ability  of  detecting  a  nonmoving  point  target  in  the  pres¬ 
ence  of  signal-dependent  reverberation  and  colored  ambient 
noise  whose  power  spectral  densities  are  known.  An  ana¬ 
lytical  solution  first  reported  by  Kooij  is  corrected  and  ex¬ 
tended.  A  simple  example  for  the  white  ambient  noise  case 
is  included  to  provide  insight  into  the  transmit  signal  opti¬ 
mization.  This  example  shows  that  the  effect  of  the  solution 
is  to  “pre-emphasize”  the  transmit  signal,  which  results  in 
the  whitening  of  the  signal-dependent  reverberation  power 
spectral  density. 

1.  INTRODUCTION 

When  an  active  sonar  or  radar  system  transmits  a  short  pulse, 
it  receives  an  echo  from  each  of  any  targets  present,  and 
from  the  many  random  scatterers  present  in  the  vicinity  of 
the  source.  The  summation  of  the  returns  from  these  ran¬ 
dom  scatterers  is  known  as  reverberation  in  the  sonar  lit¬ 
erature,  and  known  as  clutter  in  the  radar  literature.  This 
signal-dependent  reverberation  or  clutter  often  renders  the 
detection  of  any  target  echoes  present  much  more  difficult. 
This  is  especially  true  if  a  target  does  not  maintain  an  ap¬ 
preciable  radial  motion  with  respect  to  the  active  source,  be¬ 
cause  then  its  target  echo  will  not  undergo  sufficient  Doppler 
shift  that  it  can  be  readily  detected  outside  the  frequency 
band  of  the  reverberation.  Hence,  the  detection  of  echoes 
from  nonmoving  targets  in  the  presence  of  signal -dependent 
reverberation  has  always  been  considered  a  most  difficult 
detection  problem.  The  presence  of  ambient  noise  at  the 
receiver  further  complicates  the  detection  of  target  echoes. 

An  optimal  receiver  is  usually  defined  as  one  which  max¬ 
imizes  the  probability  of  detecting  a  target  in  the  presence 
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of  interference  (i.e.,  reverberation  plus  ambient  noise)  for 
a  specified  probability  of  false  alarm.  The  solution  for  the 
design  of  such  an  optimal  receiver,  given  a  statistical  char¬ 
acterization  of  the  interference,  is  available  in  a  number  of 
standard  texts  [1,  2],  It  can  be  shown  that  the  detection  per¬ 
formance  of  the  optimal  receiver  is  a  monotonic  function  of 
the  signal  to  interference  ratio  (SIR).  When  the  interference 
is  white  in  frequency,  the  SIR  depends  only  on  the  total  en¬ 
ergy  of  the  transmit  signal.  However,  when  the  interference 
is  not  white,  and  furthermore,  signal-dependent,  the  design 
of  the  transmit  signal  can  have  a  significant  impact  upon  the 
detection  performance  of  the  optimal  receiver. 

It  was  pointed  out  two  decades  ago  [3]  that  the  per¬ 
formance  of  the  optimal  receiver  in  detecting  target  echoes 
from  a  point  target  in  signal-dependent  interference  depends 
only  upon  the  channel  power  spectral  density  (PSD),  the 
ambient  noise  PSD,  the  Doppler  frequency  shift  induced  by 
the  point  target’s  motion,  and  the  power  spectrum  of  the 
transmit  signal.  When  each  of  these  factors  is  specified,  the 
solution  to  the  design  of  the  optimal  receiver  is  determined. 
Of  all  these  factors,  the  only  one  that  is  under  the  control  of 
an  active  sonar/radar  system  designer  is  the  transmit  signal. 

For  this  reason,  the  design  of  a  transmit  signal  which 
maximizes  the  performance  of  an  optimal  receiver  in  detect¬ 
ing  a  point  target  in  signal-dependent  interference  has  been 
an  active  research  topic  for  many  years  [4,  5,  6,  7,  8,  9].  In 
particular,  Kooij  [7]  has  reported  an  analytical  solution  to 
the  problem  of  transmit  signal  optimization  in  the  presence 
of  signal-dependent  reverberation  plus  ambient  noise.  How¬ 
ever,  Kooij  assumed  that  the  ambient  noise  spectral  level 
was  white  (i.e.,  constant  over  the  transmit  frequency  band). 
Also,  he  employed  a  calculus  of  variations  approach  that  op¬ 
timized  the  transmit  signal  power  spectrum  by  solving  for  a 
local  extremum,  but  he  did  not  attempt  to  show  that  his  an¬ 
alytical  solution  was  indeed  a  global  maximum. 

Work  by  the  current  authors  corrects  and  extends  Kooij ’s 
work,  and  provides  an  analytical  solution  to  the  problem  of 
optimizing  the  transmit  signal  power  spectrum  to  maximize 
the  detection  performance  of  an  optimal  receiver  in  detect¬ 
ing  target  echoes  from  a  nonmoving  point  target  in  the  pres- 
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ence  of  signal-dependent  reverberation  and  colored  ambient 
noise.  In  addition,  this  paper  interprets  the  analytical  so¬ 
lution  through  a  simple  example,  and  provides  the  results 
of  numerical  calculations  which  illustrate  the  performance 
gain  which  may  be  achieved  under  certain  conditions. 


only  the  magnitude  of  the  Fourier  transform  of  s(t)  that  af¬ 
fects  A.  Hence,  we  seek  to  maximize  A  over  all  transmit 
signal  power  spectra  Ps(f )  that  are  bandlimited  to  IF  Hz, 
and  for  which  the  total  energy  is  constrained.  We  assume 
that  the  total  transmit  signal  energy  to  be  constrained  to  E, 
or 


2.  PROBLEM  STATEMENT 

We  wish  to  detect  a  complex  signal  As(t)  embedded  in 
wide-sense  stationary  colored  and  complex  Gaussian  noise. 
The  signal  is  known,  except  for  its  amplitude  ,4,  which  is 
modeled  as  a  complex  Gaussian  random  variable  with  zero 
mean  and  variance  a\.  The  signal  is  furthermore  assumed 
to  be  bandlimited  to  IF  Hz. 

The  interference  has  the  PSD 


+w  r+W 

\S{f)?df=  /  P,(f)df  =  E , 

-w  J-w 

and  that  Ps(f)  >  0  for  |/|  <  IF. 

In  the  next  section,  the  solution  to  this  optimization  prob¬ 
lem  is  given. 

3.  THEOREM  STATEMENT  AND 
INTERPRETATION 


P(f)  =  Pc(f)PsU)  +  Pn(f), 

for  |/|  <  W,  where  Pc(f )  is  the  PSD  of  the  random  chan¬ 
nel,  Ps(f )  is  |S(/)|2,  where  S(f)  is  the  Fourier  transform 
of  s(t),  and  P„(f )  is  the  PSD  of  the  ambient  noise. 

The  optimal  receiver  can  be  shown  to  compute  the  test 
statistic 


T  = 


r  +  W 

J-W 


+w  X(f)S*(f) 


P(f ) 


df 


where  X(f)  is  the  Fourier  transform  of  x(t),  the  received 
waveform.  This  detector  is  an  asymptotic  form,  under  the 
assumption  of  stationary  processes  and  long  observation  times 
(SPLOT)  [1].  It  is  valid  if  the  observation  time  is  much 
greater  than  the  correlation  time  of  the  noise.  Under  the 
SPLOT  assumption,  the  performance  of  the  optimal  detec¬ 
tor  is 


Theorem:  Given  the  functional 


J  = 


+ll' 


P.(f) 


-W  Pc(f)P.(f)  +  Pn(f) 


df, 


where  Pc(f)  and  Pn(f )  are  known  functions  for  \f\  <  W, 
and  E  >  0  is  known,  and  where  the  selection  of  the  function 
Ps(f )  is  subject  to  the  constraint 


r+W 


Ps(f)df  =  E, 


t-w 


and  to  the  constraint  that  Ps(f )  >  0  for  |/|  <  IF,  then  the 
function  Ps(f )  which  maximizes  J  is  given  by  the  following 
expression: 


Ps(f)opt  =  max 


'  VPn(f)/\-Pn(f) 
Pc(f) 


pd  =  p1f/{1+a\ 

where  Pd  is  the  probability  of  detection,  Pp  is  the  proba¬ 
bility  of  false  alarm,  and  A  is  the  detection  index  defined 
by 


where  A  can  be  found  by  numerically  solving  the  following 
expression: 


\/Pn(f)/X-Pn(f) 
Pc(f ) 


,0  )df  =  E. 


A  =  4 


Ps(f) 

P(f ) 


df. 


Note  that  for  a  fixed  Pp,  Pd  increases  as  A  increases. 
For  the  purpose  of  simplicity,  we  assume  that  a  \  =  1,  so 
that  the  optimal  receiver  detection  performance  is  solely  de¬ 
pendent  on  the  index: 


A  = 


Ps(f) 

Pc(f)Ps(f)  +  Pn(f) 


df 


We  can  better  understand  this  theorem  by  looking  at  a 
simple  example.  Let  us  assume  that  Pn(f)  =  Xq.  Then  the 
optimal  function  Ps{f)  must  maximize 

j  =  rw _ m _ df 

'  J-W  Pc(f)Ps(f)  +  No  ^ 

subject  to  the  constraints  Ps(f)df  =  E  and  Ps(f )  > 
0  for  |/|  <  W. 

The  theorem  in  this  case  yields  the  optimal  signal  PSD 


For  best  performance,  we  would  like  to  maximize  A 
over  all  signals.  Note  that  for  the  SPLOT  condition,  it  is 


Ps(f)opt  =  max 


l  s/X o/X  -  N0  \ 

V  Pcif)  ’  J 
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where  A  is  found  by  solving 
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yi%JX-N0 

Pc(f) 


,0  \df  =  E. 


Since  (  \/Nq/\  —  No)  must  be  a  positive  constant  to  satisfy 
the  constraint,  we  have  that 


(v^VoA  -  No) 


E 

S-Zl/Pdf)df 


and  therefore 

pm  \/NoJX  —  Np  E/Pc(f) 

sWopt  PcU)  J+Zi/PC(M 


For  example,  if  Pc(f)  =  1,  then  the  optimal  Ps(f )  = 
E/2W,  which  agrees  with  a  well  known  result  [4],  Because 
this  expression  for  Ps(f )  results  from  the  assumption  that 
Pc{f)  =  1,  we  often  denote  this  expression  as  Ps(f)fiat- 
Similarly,  we  denote  as  J fiat  the  value  of  the  functional  J 
that  results  when  the  transmit  signal  power  spectrum  Ps(f ) 
is  selected  to  be  Ps(f )  fiat <  even  when  Pc(f )  is  not  equal  to 
1. 

4.  SIMPLE  EXAMPLE  WITH  NUMERICAL 
RESULTS 


f  in  Hz 


Fig.  1.  Random  channel  PSD  Pc{f)  and  ambient  noise  PSD 
Nq  for  different  values  of  the  3  dB  bandwidth  B 


Fig.  2.  Performance  gain  =  Wlog10(Jmax  /  Jfiat) 

We  assumed  that  the  PSD  of  the  random  channel  was 
Pc(f),  where 

Pc(f)  =  0.001(1  +  ((27r/)/a)2). 

Note  that  the  3  dB  bandwidth  B  of  this  channel  is  then 
a/( 2tt).  Figure  1  is  a  graph  of  this  PSD  Pc(f)  for  several 
values  of  the  3  dB  bandwidth  B. 

For  simplicity,  we  also  assumed  that  the  energy  E  =  1, 
the  signal  bandwidth  is  W  =  1000  Hz,  and  the  white  am¬ 
bient  noise  PSD  Pn(f)  =  N0  =  0.0005.  MATLAB  was 
used  to  numerically  integrate  1  /Pdf)  over  all  frequencies 
between  —W  and  +W,  so  that  Ps{f)0pt  could  be  calcu¬ 
lated.  Next,  we  calculated  Jma;e ,  the  value  of  the  functional 
J  when  Ps(f)0pt  is  substituted  for  Ps(f).  Finally,  the  de¬ 
tection  performance  gain  was  calculated  as 

GAIN  =  Wlogio(Jmax/Jflat)- 

Figure  2  is  a  graph  of  this  detection  performance  gain 
in  dB  that  results  from  selecting  the  optimal  transmit  sig¬ 
nal  power  spectrum,  for  several  different  values  of  the  3  dB 
channel  bandwidth  B.  It  is  clear  that  significant  gains  may 
be  possible  under  certain  channel  conditions. 

Figure  3  is  a  graph  of  the  reverberation  PSD  (Pc(f)Ps(f)) 
when  Ps(f)  =  Ps{f) fiat  and  when  Ps(f)  =  Ps(f)opt- 
This  figure  shows  that  the  effect  of  using  Ps(f)opt  is  to 
whiten  the  reverberation  PSD. 

5.  CONCLUSIONS 


To  illustrate  that  the  selection  of  an  optimal  transmit  sig¬ 
nal  power  spectrum  can  improve  the  detection  performance 
of  an  optimal  receiver,  we  used  MATLAB  to  calculate  the 
gain  in  the  functional  J  for  Ps(f)  =  Ps{f)opt  over  the 
choice  of  Ps(f)  =  Ps(f)fiai- 


In  this  paper,  we  have  presented  a  new  result  for  transmit 
signal  optimization  for  active  sonar/radar.  The  result  pre¬ 
sented  corrects  and  extends  an  analytical  solution  originally 
reported  by  Kooij.  This  result  is  only  valid  for  detecting 
target  echoes  from  nonmoving  point  targets  embedded  in 
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Fig.  3.  Ambient  noise  PSD  N(h  and  reverb  spectrum 
Pc(f)PsU)  when  Ps(f)  =  Ps{f) flat  and  when  Ps(f)  = 
Ps(f)optJoiB  =  10 

signal-dependent  reverberation  and  colored  ambient  noise. 
A  simple  example  was  included,  to  demonstrate  that  the  ef¬ 
fect  of  the  transmit  signal  optimization  can  lead  to  signifi¬ 
cant  detection  gains.  Also,  it  was  noted  that  for  white  am¬ 
bient  noise,  the  optimal  transmit  signal  is  chosen  so  that 
the  resulting  signal-dependent  reverberation  will  have  a  flat 
PSD  over  the  transmit  frequency  band. 
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